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The dynamical behavior of a higher-order cubic Ginzburg-Landau equation is found to include 
a wide range of scenarios due to the interplay of higher-order physically relevant terms. We find 
that the competition between the third-order dispersion and stimulated Raman scattering effects, 
gives rise to rich dynamics: this extends from Poincare-Bendixson-type scenarios, in the sense that 
bounded solutions may converge either to distinct equilibria via orbital connections, or space-time 
periodic solutions, to the emergence of almost periodic and chaotic behavior. One of our main 
results is that the third-order dispersion has a dominant role in the development of such complex 
dynamics, since it can be chiefly responsible (i.e., even in the absence of the other higher-order effects) 
for the existence of the periodic, quasi-periodic and chaotic spatiotemporal structures. Suitable low¬ 
dimensional phase space diagnostics are devised and used to illustrate the different possibilities and 
identify their respective parametric intervals over multiple parameters of the model. 

PACS numbers: 02.30.Jr, 05.45.-a 


I. INTRODUCTION 

Nonlinear evolution equations are often associated 
with the theory of solitons and integrable systems [T]. A 
prime example is the nonlinear Schrodinger (NLS) equa¬ 
tion which constitutes one of the universal nonlinear evo¬ 
lution equations, with applications ranging from deep wa¬ 
ter waves to optics [5]. Remarkable phenomena are also 
exhibited by its higher-order variants, emerging in a di¬ 
verse spectrum of applications, such as nonlinear optics 
[5] , nonlinear metamaterials [1] , and water waves in finite 
depth On the other hand, dissipative variants of 

NLS models incorporatiirg gain and loss have also been 
used in optics [5], e.g., in the physics of mode-locked 
lasers pun] (see also the relevant works mmn) and 
polariton superfluids [T3| - see, e.g.. Ref. m for various 
applications. Note that such dissipative NLS models can 
be viewed as variants of the complex Ginzburg-Landau 
(CGL) equation, which has been extensively studied, es¬ 
pecially in the context of pattern formation in far-from- 
equilibrium systems Ha- 

Dissipative nonlinear evolution equations (incorporat¬ 
ing gain, loss, external driving, or combinations thereof) 
may exhibit (and potentially be attracted to) low¬ 
dimensional dynamical features, such as: (a) one or more 
equilibria (and orbits connecting them), (b) periodic 
orbits, (c) quasi-periodic orbits or (d) low-dimensional 
chaotic dynamics m- The availability of the dynamical 
scenarios (a)-(d) depends on the effective dimensional¬ 
ity of the low dimensional behavior; one-dimensionality 
only allows fixed points, planar systems governed by the 
Poincare-Bendixson (PB) theorem [1^1 can also feature 


periodic orbits, while higher dimensions allow for quasi- 
periodic or chaotic dynamics. Various prototypical par¬ 
tial differential equation models have demonstrated a 
PB-type behavior as an intermediate bifurcation stage 
in the route to spatiotemporal chaos. Examples include 
the Kuramoto-Sivashinsky m and complex Ginzburg- 
Landau (CGL) equations; regarding the CGL model, 
which is of primary interest in this work, we refer to 
the seminal works |18] for the spatiotemporal transition 
to chaos. In addition to the above autonomous systems, 
spatiotemporal chaos was also found in non-autonomous 
ones, due to the interplay between loss and external 
forces, such as the damped-driven NLS [TOMT] (where 
the hyperbolic structure of the underlying integrable NLS 
is a prerequisite [55]) and the sine-Gordon [^ system. 

In this work, we focus on the role of higher order ef¬ 
fects, and investigate the possibility of bifurcation phe¬ 
nomena leading to the existence of the above proto¬ 
typical examples of low-dimensional dynamics in an au¬ 
tonomous, physically important higher-order CGL-type 
model. This model, is motivated by the higher-order 
NLS equation that is commonly used, e.g., in studies of 
ultrashort pulses in optical fibers |3], but also incorpo¬ 
rates (linear or nonlinear) gain and loss; it is, thus, a 
physically relevant variant of a higher-order cubic CGL 
equation - without the diffusion term. Note that higher- 
order versions of the CGL equation, have only recently 
started attracting attention [57], while extended second- 
order CGL models have been extensively studied in var¬ 
ious contexts previously 0 m [15]. In particular, we 
refer to the pioneering work [lOj . followed by the impor¬ 
tant contributions [HI Ha, which revealed the existence 
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of the aforementioned low-dimensional dynamical scenar¬ 
ios for second-order quintic CGL models. The results of 
[inHi2] were established by numerical and even analyt¬ 
ical reductions to suitable finite dimensional dynamical 
systems, capturing the long-time dynamics of the original 
infinite dimensional one. Notably, the revealed dynami¬ 
cal scenarios were associated with a variety of novel local¬ 
ized structures (known as pulsating solitons). However, 
a crucial feature of these models as acknowledged, e.g., 
by the authors of [10] was the presence of a higher order 
(quintic) nonlinearity. This is a key trait distinguishing 
that set of works from the present one where only cubic 
nonlinearity is employed, yet the presence of higher or¬ 
der effects, most notably third order dispersion as we will 
see below, plays a catalytic role in the emergence of the 
relevant phenomenology. 

More specifically, we should point out that the third- 
order cubic CGL model that we consider herein is essen¬ 
tially different from the second-order cubic-quintic model 
discussed in [Mil], not only from a mathematical but 
also from a physical point of view: indeed, in the context 
of optics, the model considered in the latter works refers 
to propagation of short pulses, in the picosecond regime, 
in media featuring saturation of the nonlinear refractive 
index, while the model we consider here is relevant to 
propagation of ultrashort pulses in the sub-picosecond or 
femtosecond regimes [3]. For this reason, our model in¬ 
cludes third-order dispersion and higher-order nonlinear 
effects, that appear naturally as higher-order corrections 
of the usual NLS model in the framework of the reductive 
perturbation method. In that regard, and if gain and loss 
are also incorporated, it is important to ask if, and how, 
the physically important (in the femtosecond time-scale) 
higher-order effects may be responsible for tracing a path 
to complex dynamics, as a result of the potential break¬ 
ing of the homoclinic structure of the unperturbed NLS 
counterpart. 

The main findings of our investigations are the follow¬ 
ing. First, we show that the incorporation of the gain 
and loss terms gives rise to the existence of an attrac¬ 
tor; a rigorous proof is provided, based on the interpre¬ 
tation of the energy balance equation and properties of 
the functional (phase) space on which the problem de¬ 
fines an infinite-dimensional flow. The structure of the 
attractor is then investigated numerically. Given that our 
model is characterized by six free parameters (which ren¬ 
ders a systematic investigation of their role a nontrivial 
task), we opt to keep four parameters fixed, with values 
suggested by the physics of ultrashort optical pulses [3], 
and vary the remaining two. In particular, we vary the 
coefficients of the third-order dispersion and the higher- 
order nonlinear dissipation, accounting for the stimulated 
Raman scattering (SRS) effect (more important reasons 
for this choice will become apparent below). We find 
that, for sufficiently small SRS coefficient, variations of 
the third-order dispersion strength give rise to a tran¬ 
sition path from dynamics reminiscent of PB, including 
orbital connections between steady states of high multi¬ 


plicity and convergence to limit cycles, to invariant tori 
or even chaotic attractors. However, when the SRS ef¬ 
fect becomes stronger, the above scenarios are screened 
by convergence to steady-states. It is highlighted that 
the third-order dispersion is found to be chiefly respon¬ 
sible for a dynamical transition from periodic, to quasi- 
periodic and eventually to chaotic structures. Therefore, 
our results show that higher (third)-order dispersion and 
dissipative (SRS) effects are important mechanisms for 
the emergence of complex spatiotemporal transitions in 
GGL models. 

Our presentation is organized as follows. In Section H, 
we present the model, and discuss the existence of a limit 
set (attractor). Details on the proof of such a limit set 
are given in the Appendix]^ The structure of the attrac¬ 
tor is then investigated numerically in Section HI. We 
thus reveal the emergence of all dynamical scenarios and 
corresponding regimes of complex asymptotic behavior. 
Finally, Section IV summarizes our findings. 


II. MOTIVATION AND PRESENTATION OF 
THE MODEL 

Our model is motivated by the following higher-order 
NLS equation: 

idtu — -\- \u\‘^u = -\- iqLdx{\u\‘^u) 

+ {crR + w)udx{\u\'^), ( 1 ) 

where u{x, t) is a complex field, subscripts denote partial 
differentiation, j5, /x, v and ur are positive constants, 
while s = ±I denotes normal (anomalous) group veloc¬ 
ity dispersion. Note that Eq. Q can be viewed as a 
perturbed NLS equation, with the perturbation (in case 
of small values of relevant coefficients) appearing in the 
right-hand side (see, e.g.. Refs. [3] and discussion below). 

Variants of Eq. Q appear in a variety of physical con¬ 
texts, where they are derived at higher-order approxima¬ 
tions of perturbation theory [the lowest-order nonlinear 
model is simply the NLS equation in the left-hand side of 
Eq. Q]. The most prominent example is probably that 
of nonlinear optics [3] . In this case, t and x denote prop¬ 
agation distance and retarded time, respectively, while 
u{x,t) is the electric field envelope. While the unper¬ 
turbed NLS equation is sufficient to describe optical pulse 
propagation, for ultra-short pulses third-order dispersion 
and self-steepening (characterized by coefficients /3, /x and 
R, respectively) become important and have to be in¬ 
corporated in the model. Similar situations also occur 
in other contexts and, thus, corresponding versions of 
Eq. (§ have been derived and used, e.g., in nonlinear 
metamaterials [3] , but also in water waves in finite depth 
Moreover, in the context of optics, and for rela¬ 
tively long propagation distances, higher-order nonlinear 
dissipative effects, such as the SRS effect, of strength 
CTfl > 0, are also important [5]. 
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In addition to the above mentioned effects, our aim is 
to investigate the dynamics in the framework of Eq. Q , 
but also incorporating linear or nonlinear gain and loss. 
This way, in what follows, we are going to analyze the 
following model: 

idtu —-d^u + \u\'^u = iju + i5\u\'^u + ifidx{\u\'^u) 

+ il3dlu+{aR + iiy)udxi\u\‘^), ( 2 ) 

which includes linear loss (7 < 0 ) [or gain (7 > 0)]. 
These effects are physically relevant in nonlinear optics 
[31IH1IT1]: indeed, nonlinear gain (6 > 0) [or loss {S < 0)] 
may be used to counterbalance the effects from the lin¬ 
ear loss/gain mechanisms and can potentially stabilize 
optical solitons - see, e.g., Refs. [Ml US]- As is also ex¬ 
plained below, here we focus on the case of linear gain, 
7 > 0 , and nonlinear loss, 5 < 0 , corresponding to a 
constant gain distribution, and the intensity-dependent 
two-photon absorption, respectively (see, e.g.. Refs. 

Obviously, the presence of gain/loss renders Eq. Q 
a higher-order cubic CGL equation (cf. recent studies 
m on such models). Note that in Eq. ([^, diffusion is 
absent: such a linear term would be of the form iDd^u 
(D =const.), and would appear in the right-hand side of 
Eq. (§ to account for the presence of spectral filtering 
or linear parabolic gain {D > 0) or loss {D < 0) (TUl- 
[H]. Instead, the equation only features linear dispersion 
through the term proportional to s in the left hand side 
of Eq. ([ 2 ]) . 

The gain/loss effects are pivotal for the dissipative na¬ 
ture of the infinite-dimensional flow that will be defined 
below. This dissipative nature is reflected in the exis¬ 
tence of an attractor, capturing its long-time dynamics; 
nevertheless, as we will show below, the structure of the 
attractor is determined by the remaining higher-order ef¬ 
fects. 

Here, we focus on the case s = 1, and supplement 
Eq. (U) with periodic boundary conditions for u and its 
spatial derivatives up to the-second order, namely: 

u{x + 2L,t) = u{x,t), and 
di{x + 2L,t) = di{x,t), j = l,2, 

V {x,t) £ R X [0,r], for some T > 0, where L > 0 is 
given. The initial condition 

u{x,0) = uo(x), Va:GM, (4) 

also satisfies the periodicity conditions Here, we 
should mention that the periodic boundary conditions 
that we consider here are also motivated by the context of 
optics. Recalling that the roles of space and time are in¬ 
terchanged in the latter context, we note that, indeed, in 
optical cavities (e.g., ones for lasers), the period L would 
account for the (retarded) time it takes light to traverse 
to the laser cavity once and, thus, the boundaries rep¬ 
resent the same point in the real space-time (see, e.g.. 
Refs. [MH3I])- In this context, the dynamics that will be 
analyzed below are relevant to the dynamical transitions 


and the observation of chaotic optical waveforms in fiber 
ring lasers [3fi] . 

As shown in Ref [35]) all possible regimes except 7 > 0, 
(5 < 0 , are associated with finite-time collapse or decay. 
Eurthermore, a critical value 7 * can be identified in the 
regime 7 < 0 , (5 > 0 , which separates finite-time collapse 
from the decay of solutions. On the other hand, for 7 > 0, 
(5 < 0, we prove in Appendix [^ the existence of a limit 
set (attractor) uj^Uq), attracting all bounded orbits ini¬ 
tiating from arbitrary, appropriately smooth initial data 
uq (considered as elements of a suitable Sobolev space). 
In the next Section, we will show numerically that the 
attractor uj{uo) captures the full route from PB-type dy¬ 
namics to quasi-periodic or chaotic dynamics. 


III. NUMERICAL RESULTS 

The structure of the limit set uj(uo), is investigated 
by numerical integration via a high-accuracy pseudo- 
spectral method. In our simulations, we fix the half 
length of fl to L = 50, and the ratio —j/S to be of 
the order of unity, and thus fix 7 = 1.5 and 5 = — 1. 
This choice, which stems from the fact that this ratio de¬ 
termines the constant density steady-state (see below), 
will be particularly convenient for illustration purposes. 
Furthermore, motivated by the fact that, in the context 
of optics, parameters describing the higher-order effects 
take, typically, small values |3|, we fix /x = = 0.01, 

while third-order dispersion and SRS strengths, /3 > 0, 
(Tr > 0, are varied in the intervals [0, 1] and [0, 0.3], 
respectively. 

Obviously, the above choice is merely a low¬ 
dimensional projection of the full 6 -dimensional param¬ 
eter space. Nevertheless, since our scope here is to illus¬ 
trate the role of higher-order effects on the emergence of 
complex dynamics in Eq. ([^, we will show below that 
the variations of /3 and aR alone do offer a clear physical 
picture in that regard. To be more specific, the choice 
of those particular parameters stems from the follow¬ 
ing facts. First, third-order dispersion is the sole linear 
higher-order effect, which is important also in the linear 
regime (as it modifies the linear dispersion relation). Sec¬ 
ond, the stimulated Raman scattering effect is the first 
higher-order dissipative effect and, as such, is expected to 
play dominant role in the long-time nonlinear dynamics 
of the system. 

Naturally, the nontrivial task (as also highlighted 
above) of investigating the full parameter space is in¬ 
teresting and relevant in its own right, yet it is beyond 
the scope of this work. 

In our simulations, the limit set oj(uq) will be vi¬ 
sualized by projections of the flow to suitable 2D or 
3D spaces, defined by 7^2 = {(A, F) € M^}, and V 3 = 
{(X,Y,Z) €R^}. Here, X(t) = \u{xi,t)\^, Y{t) = 
|u(x 2 ,t)|^, Z{t) = \u{x 3 ,t)\'^, for arbitrary spatial coor¬ 
dinates Xi, X2,X3 G fl. 
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Figure 1: (Color Online) The scenario aj(uo) = {(fib}- Left 
panel: convergence to the fixed point A. Right Panel: the 
fixed point A as a limit circle of radius 


A. Steady-state and orbital connections regime 


First, we use cw initial data, 


Uq{x) = e exp 


Kttx\ 


e4>K 


of amplitude e > 0 and wave-number K > 0, which is an 
element of the ID-linear subspace 

Vk = {u £ L^{n) : u = e(/)if(x), e > 0} 

of L^(0). Here we should note that there exists a cw 
state which is an exact solution of Eq. § ; this solution 
is generically subject to modulational instability (MI) 
[85] (so-called Benjamin-Feir instability in the context 
of deep water waves |3B] ). The exact cw solution, as well 
the relevant MI analysis are presented in Appendix B. 
However, such analysis is not capable of providing any 
insights on the long-time dynamics of the solutions. In¬ 
deed, although it can be used as a means to understand 
the destabilization of the cw steady-state, it does not of¬ 
fer any information regarding the long-time behavior and 
the states the system passes through. As we show below, 
the intricate dynamics that emerge, cannot be fully un¬ 
derstood in the framework of the MI picture. 

Using the above cw initial data, and varying an > 0, 
we find that uj^uq) is an equilibrium state. Specifically, 
there exists a critical wave number iFmax such that: for 
K < ATmax, w(uo) = (f>b, Ls., a steady-state of constant 
density for K > u}{uo) = $p, i.e., a 

steady-state of spatially periodic density. We find that 
ATmax decreases as an increases: ii an = 0,0.1, 0.2,0.3, 
and fi = 0.02, then ATmax = 16,13,10, 5, respectively. 

The dynamical scenario a;(uo) = {4>b} for /3 = 0.02, 
an = 0.3 and AT = 4 is illustrated in Fig. [T] The projec¬ 
tion of the cw equilibrium (f>b to the 2D space U 2 is the 
fixed point A = (|^f,|2, |(/,b|2) = = (1.5,1.5). 

The right panel of Fig. illustrates the convergence of 
the projected linear orbits to A, associated to the choice 
of spatial coordinates Xi = 5 , 0:2 = 10. The dashed blue 
(continuous red) line is the projection of the flow for the 



Figure 2: (Color Online) The scenario cj(uo) ~ {“Lp}. Upper 
panels: density snapshots at times (a) t « 500, (b) t ~ 683, 
(c) t « 700. Bottom panels: orbital connections O A —>■ B 
in 2D (left) and 3D (right) spaces. 


cw with e = 3 (e = 0.01); the arrows indicate the direc¬ 
tion of the 2D-projection of the flow. The cw steady state 
(f>b is an element of Vk , and only differs in amplitude from 
the initial condition. Hence, Vk defines a stable linear 
subspace for A. The right panel of Fig. [^visualizes the 
steady state (j)b as a limit circle A of radius = v^LS, 
in the 2D space (Re(M(0, t)), lm(u(0, t)). The limit circle 
corresponds to the rotating linear oscillations of the real 
and imaginary parts of the solution u. Effectively in this 
case, the solution preserves its plane-wave form but its 
amplitude, say h{t), satisfies a Bernoulli equation. This 
can be seen as follows: we substitute u{x,t) = W{t)e'^^ 
in Eq. ([^ and obtain the following equation for the time- 
dependent amplitude W{t) (as usual overhead dots de¬ 
note time derivatives): 

iW - + \W\‘^W = i'yW + iSlWl^W 

+ PK^W - ^iKW. ( 5 ) 

Then, taking W{t) = h(t) exp[i(Ar^/2 —/3Ar^-|-^A')t], we 
obtain from Eq. the Bernoulli equation: 

h = ^h + 5hP. 


Thus, for /i(0) = e, 

lim h?{t) = -% = \(t)b?- 

t—>-oo 0 


Next, consider the scenario a;(uo) = {4>p}, for /3 = 0.02 
an = 0.3, and AT = 5, illustrated in Fig. [^ The upper 
panel shows density snapshots for a cw-initial condition 
with e = 0.01. The solution has reached the cw-steady 
state (j)b exponentially fast, but at t « 500 the instability 
of the state (j)b emerges. Although transient oscillations 
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Figure 3: (Color Online) The dynamics scenario tu(uo) = L, 
i.e., a space-time periodic traveling wave. Upper panels: den¬ 
sity snapshots at times (a) t « 135, (b) t « 150, (c) t « 180. 
Bottom panels: convergence O —>■ A —> L, the limit cycle in 
2D (left) and 3D (right) spaces. 

of increasing amplitude occur (cf. snapshot at t = 683) 
due to the linear gain 7 > 0, the nonlinear loss 5 < 
0 prevents collapse of the solution. After t « 685, we 
observe convergence to the new steady state <l>p (reached 
at < « 700), whose profile remains unchanged till the end 
of integration {t = 3000). The orbital connection, via 
the transient dynamics, between steady states (j)b and 
is illustrated in the projections of the flow on the 
spaces 7^2 and Vs - cf. bottom left and right panels 
of Fig. respectively, for xi — 0 and X 2 = 4.5. In 
2D, B « (1.5,0.15) is the new fixed point, while in 3D, 
A = (1.5,1.5,1.5) and B « (1.5, 0.15,1.16). The infinite¬ 
dimensional orbital connection: 

{ 0 } (unstable) {(j)^} (unstable) {d’p} = w(uo), 

where Oi and O 2 denote the orbits connecting the steady 
states, is projected to the 2D and 3D-connections: 

O (unstable) {A} (unstable) {B}. 

The projected orbits highlight the spiraling of the stable 
manifold of the limit point B around the unstable linear 
subspace of O = (0,0, 0) connecting O and A. The con¬ 
nection was found to be stable with respect to variations 
of e - cf. linear dashed blue (continuous red) converg¬ 
ing orbit in the bottom left panel, corresponding to a 
cw-initial condition of amplitude e = 2 (e = 0 . 01 ). 

B. Space-time periodic (limit-cycle) regime 

Increasing /3, for = 0.01, we observe the birth of yet 
another feature, namely traveling space-time oscillations. 
The upper panel of Fig. shows density snapshots, for 


Figure 4: (Color Online) The dynamics scenario tu(uo) = L, 
in the presence of third-order dispersion only, namely, for /3 = 
0.02 and ^ = v = an = Q. Upper panel: density snapshots 
at times (a) t « 400, (b) t ~ 420 and (c) t ~ 450, for a 
single-cw initial condition of A = 15, e = 0.01, and /3 = 0.02. 
Bottom left and right panels: convergence O — > A — )■ L, i.e., 
the stable limit cycle, in the 2D and 3D spaces. 


a cw initial condition of A = 5, e = 0.01, and /? = 0.55. 
Now, instability of the steady-state (/>(,, leads to the birth 
of a stable, traveling space-time periodic solution, whose 
profile is shown for t = 180 (arrow indicates propaga¬ 
tion direction). The projections, for xi = 0, X 2 = 5 
and 2:3 = 10 , on V 2 (bottom left panel) and V 3 (bottom 
right panel), visualize the periodic solution as a limit cy¬ 
cle L, i.e., a periodic orbit. The continuous blue (dashed 
red) linear orbit shown in the bottom left panel corre¬ 
sponds to the cw-initial condition of A = 5 and e = 3 
(e = 0.01), highlighting the stability (i.e., attracting na¬ 
ture) of the limit cycle with respect to e. Specifically, 
for fixed an = 0.01 and A > 4, there exists an interval 
Ip,K = [/3min(Ar),/3inax(A:)], such that for some /3 G Ip^K, 
the initial condition may converge to a space-time peri¬ 
odic solution; e.g., for A = 5, « [0.5,0.57], while 

for A = 20, 7 / 3,20 ~ [0.7,1.2]. On the other hand, when 
/3 i the initial condition converges to a steady state. 
Evidently, the structure of the limit set w(uo) for Eq. (|^, 
consisting either of distinct equilibria and orbits connect¬ 
ing them, or of a limit cycle, is reminiscent of scenarios 
associated with PB dynamics. 

It is important to remark that third-order dispersion 
plays a critical role in this scenario of uj(uo) = L, as 
it can be solely responsible for the emergence of a limit 
cycle. Indeed, Fig. shows the dynamics for a cw-initial 
condition of A = 15 and amplitudes as in Fig. but for 
/3 = 0.02 and an = yL = v = 0. Furthermore, the third- 
order dispersion alone, can also give rise to even more 
complex behavior (see below). 
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Figure 5: (Color Online) Birth of a chaotic attractor cj(uo) = 
S. Transition from the instability of the cw-steady state A, 
to quasiperiodic, and to chaotic behavior for t € [0,330]. 


C. Quasi-periodic and chaotic regime 


The interval Ip^K rnay be partitioned to sub-intervals 
where quasi-periodic, or even chaotic behavior emerges. 
Figure shows the 3D-projection of the flow on for 
xi = 5, X 2 = 10, xz = 15, t € [0,350], /3 = 0.52, 
CTfi = /r = r/ = 0.01, for a cw of e = 0.01 and K = 5. We 
observe the birth of quasi-periodic orbits from the insta¬ 
bility of the steady-state and the transition to chaotic 
behavior manifested by their trapping to a chaotic attrac¬ 
tor S. 

The upper left panel of Fig. shows part of a chaotic 
orbit in S, for t S [180,200], and /3 = 0.5 « /3min(5). 
The first two snapshots of the bottom panel show profiles 
of the solution corresponding to points PI and P2, for 
t = 150 and t = 165. The “windings” of the chaotic or¬ 
bits are evident in the upper left panel of Fig.[^ similarly 
also to the bottom right panel of Fig.j^ The chaotic be¬ 
havior manifests itself in the time-fluctuating amplitude, 
the changes in the waveform’s spatial periodicity, and in 
the propagation direction of the chaotic traveling wave. 

The interval I/^^k = [l3rain{K), can be 

partitioned in the following sub-intervals: a chaotic 
d/3, if, c = [/3min(dd),/?ch(dd)], a quasi-periodic Ip,K,q = 
{/3ch{K),/3ic{K)), and a limit-cycle one = 

[/3ic(iF),/3max(dd)]- Let /3inin(dd) be the critical value for 
the onset of the quasiperiodic behavior and the tran¬ 
sition to the chaotic regime. Then, as f3 ^ PchiK), 
the chaotic features are less evident and emerge at later 
times. Chaotic orbits still exist for /3 = P^hiK). For 
P > Pc\i(K), solutions remain quasi-periodic, and the or- 



Figure 6: (Color Online) Top left panel: a chaotic path in S 
for t £ [180, 200]. Top right panel: projection in 3D-space Vs 
of the invariant torus-like set Q for t £ [1800,2000]. Bottom 
panels: chaotic waveforms, corresponding to points PI at 
time t ~ 150 (left) and P2 at time t ~ 165 (middle), and a 
quasi-periodic solution in Q at time t ~ 1900 (right). 


bit is trapped within an invariant torus-like set Q. For 
K = 5, we find that /3ch(5) « 0.53. The projection on V 3 
of Q for /3 = 0.54 > /3ch(5), is shown in the upper right 
panel of Fig. The orbit is plotted for t G [1800, 2000], 
and the profile of a quasi-periodic solution within Q at 
t = 1900 is shown in the third snapshot of the bottom 
panel. The set Q persists as long as P < Pci{K). When 
Pic{K) < P < /3max(dF), the set Q is replaced by a limit 
cycle. For AT = 5, we find the following sub-intervals 
of d/ 3,5 ~ [0.5,0.57]: the chaotic 1 / 3 ^ 5 ,c ~ [0.5,0.53], the 
quasi-periodic Ii 3 , 5 ,q ~ (0.53,0.55), and the limit-cycle 
d/ 3 , 5 , /c ~ [0.55,0.57]. For AT = 5, the above sub-intervals 
were detected with accuracy 10“^: for P = 0.549, the set 
Q persists, while for P = 0.55, the initial state is trapped 
on the limit cycle. 


D. Numerical bifurcation diagrams 

The richness of the dynamics can be summarized in 
a bifurcation diagnostic (“diagnostic I”), namely the 
||w||oo bifurcation diagram, shown in the upper panel 
of Fig. The bifurcation curve [continuous (blue) line] 
illustrates the variations of the [|m[ [co-norm of the solu¬ 
tions, defined as 

ll^llcJO 01 5 dl [ A, A] X [0,A,^ax]; 

where T^ax denotes the end of the interval of numeri¬ 
cal integration [0,Tmax]; the third-order dispersion coef¬ 
ficient is P G [0,1], while the rest of parameters are fixed 
to the values an = fi = v = 0.01, and for the cw-initial 
condition we use e = 0.01 and AT = 5. The system was 
integrated until r^ax = 3000. The branches AB and FG 
correspond to the intervals P G [0, 0.18) and P G (0.57,1] 
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respectively, and are associated with the dynamical sce¬ 
nario a;(uo) = {4>b}, i-e., the convergence to the steady- 
state of constant density, |^f,p = — The intersection of 
the bifurcation curve with the auxiliary “separatrix” B, 
at /3 « 0.18, designates the transition to the equilibrium 
metastability region BC [light grey (pale yellow) shaded 
area], in the interval /? S [0.18,0.5). The fluctuations of 
the bifurcation curve are associated with metastable dy¬ 
namical scenarios between distinct states. One such sce¬ 
nario may refer to the orbital connections between steady 
states mentioned above; another one, may correspond to 
a transition from unstable periodic orbits to chaotic os¬ 
cillations, and an eventual convergence to a steady state. 
These scenarios are followed by drastically different tran¬ 
sient dynamics characterizing these connections. 

As a first example, we note the metastable transition 
- at 13 = 0.3 [vertical dashed (red) line] - between three 
distinct steady-states El —>■ E2 E3 (with El marking 
the steady-state of constant density, |(/>b|^ = — |). The 
third row panels of Fig. show density profiles of these 
steady states. A second example, refers to the transition 
from an unstable periodic orbit PO (which emerges from 
the instability of the steady-state (j)b), to chaotic oscilla¬ 
tions CO and the convergence to the final steady-state 
E3; this transition occurs for /? = 0.47 [horizontal dashed 
(black) line]. Density profiles during this transition are 
shown in the fourth row panels of Fig.[^ For the first ex¬ 
ample, the ultimate state E3 is reached at t « 103, and 
the solution remains unchanged until the end of integra¬ 
tion, while for the second example, the ultimate state E3 
is reached at t « 217. 

The intersection of the bifurcation curve with the sec¬ 
ond auxiliary separatrix C, with an almost vertical slope, 
is associated with the transition to the chaotic region 
CD (grey-shaded area), corresponding to the interval 
I/ 3 , 5 ,c ~ [0.5,0.53]. The sudden jump of the bifurca¬ 
tion curve (with an infinite slope) at the intersection 
with the separatrix D designates the entrance into the 
quasi-periodic regime DE [dark (pale red) shaded area], 
associated with the interval // 3 , 5 ,g « (0.53,0.55). This 
region is magnified in the second row panel of Fig.[^ On 
the other hand, the next steep jump at the intersection 
with the separatrix E (also magnified in the second row 
panel of Fig. depicts the entrance to the space-time 
periodic regime EF [grey (pale green) shaded area], asso¬ 
ciated with the limit-cycle interval 7 / 3 , 5 ,/c ~ [0.55,0.57]. 
The limit-cycle branch bifurcates from the intersection 
with the separatrix F beyond which the branch of the 
constant-density steady-state FG is traced. Another bi¬ 
furcation diagnostic (“diagnostic II”) that we use herein, 
is the one associated with the variation of the quantity 

l|u(7’max)|la = ^ y \u{x,T^^^)\'^dx 


with respect to (3. For sufficiently large T^ax, | |' »(Tm ax)||a 
could be thought of as the superior limit of Eq. (A2). The 
drawback in the above diagnostic is that the transient 
dynamics are hidden (for sufficiently large Tmax)] more 





Figure 7: (Color Online) Top panel: j3 — ||u||oo bifurcation 
diagram (Diagnostic I), for fixed (jr — ^ = u = 0.01, and 
the cw-initial condition of e = 0.01 and A = 5. Second row 
panel: Magnification of the quasi-periodic region DE shown 
in the top panel. Third row panels: Profiles of the distinct 
steady states involved in the orbital connection El —>■ E2 —>■ 
E3 occurring at /? = 0.3 in the metastable region BC. The 
system is at rest in the steady-state El for 5 < t < 35, in 
E2 for 60 < t < 68, and in E3 for 100 < t < 3000-the 
end of integration. Fourth row panels: Transition from an 
unstable periodic orbit PO to chaotic oscillations CH which 
are eventually damped to the steady state E3. The unstable 
periodic orbit PO survives for 62 < t < 105, and the chaotic 
orbit for 120 < t < 203. The system is at rest in the steady 
state E3 for 217 < t < 3000-the end of integration. 






































Figure 8: (Color Online) Top panel: j3 — ||w(Tnax)||a bifur¬ 
cation diagram (Diagnostic II), for fixed an = = v — 0.01, 

and the cw-initial condition of e = 0.01 and K = 5. 


generally, the result strongly hinges on the selection of 
Tmax, but not necessarily strongly on the evolution for 
earlier or mirroring that for later times. Nevertheless, 
for sufficiently large T^ax, it can be particularly useful 
in detecting convergence to different steady-states, e.g., 
uj{uo) = {4>b} or a;(uo) = {‘hp}, via metastability. Fur¬ 
thermore, it is also able to detect regimes of more com¬ 
plex behavior, similarly to the ||m||oo- diagnostic. Fig- 
urej^shows the/3—||u(Tniax)||a bifurcation curve [contin¬ 
uous (red) line], for Tmax = 3000; the rest of parameters 
are as in Fig. The four shaded regions correspond to 
the same distinct dynamical regimes that were detected 
in the /3 — ||u||oo bifurcation diagram of Fig. The 
horizontal straight lines 

||^x(Tmax)||^ = 1.5=-^ 

in the regions AB and FG show that, in these regimes of 
/3, solutions converge to the steady-state 4>b- The inter¬ 
section of the bifurcation curve with the auxiliary “sep- 
aratrix” B, at /3 « 0.18, still designates the transition to 
the equilibrium metastability region BC. However, the 
new horizontal straight line ||u(Tmax)||a = 0.68 clearly 
shows that, after the transient metastability dynamics, 
the solution favors a particular steady-state of conver¬ 
gence, namely E3 for these parameters. 

It is now useful to compare Diagnostics I and II. First 
we note that the comparison between the two in the 
metastability regime BC, reveals that far-from-equilibria 
transient dynamics are only identihed by the fluctuations 
in the f3 — ||m||oo curve (Diagnostic I) - and not in the 
/3 — ||u(Tmax)|la (Diagnostic II). These fluctuations can 
be understood by the fact that ||u||oo may be reached 
at a certain instant, S [0,Tmax] and also by noting 
that, in general, ||u||oo max_i<a;<i|<I>(a:)| [i.e., the 

||u||oo-norm of a steady-state $(a:)]. Diagnostic II, on 
the other hand, reveals that in the metastability regime 



Figure 9: (Color Online) Top panel: cr/j —||w(Tniax)||a bifurca¬ 
tion diagram (Diagnostic II), for fixed /3 = 0.52, ^ = n — 0.01, 
and the cw-initial condition of e = 0.01 and K = 5. Bottom 
left panel: A chaotic path for t G [600,650], when (3 = 0.53, 
Or = h = n — 0 and the initial condition is as in the top 
panel. Bottom right panel: chaotic waveforms corresponding 
to points PI at t « 600 (top), P2 at t « 625 (middle), and 
P3 at t ~ 650. 


BC, the dynamics favors a distinct steady-state (as men¬ 
tioned above) - a fact that cannot be captured by Diag¬ 
nostic I. As far as the other regimes are concerned. Di¬ 
agnostic II can also capture the transition to the chaotic 
regime CD, indicated by the intersection of the bifurca¬ 
tion curve with the auxiliary separatrix C, as well as by 
its large rapid fluctuations within region CD. The sud¬ 
den jump of the bifurcation curve at the intersection with 
the separatrix D designates the entrance into the quasi- 
periodic regime, portrayed by the small, almost horizon¬ 
tal branch of quasi-periodic solutions within region DE. 
Note that the transition to the quasi-periodic regime is 
much more apparent in the Diagnostic II than in Diag¬ 
nostic I. The intersection of the bifurcation curve with 
the separatrix E (at a point where the curve has a lo¬ 
cal minimum in the region DF), is again associated with 
the entrance to the space-time periodic regime EF (corre¬ 
sponding to the branch of space-time periodic solutions). 
This branch bifurcates from the straight line FG (perti- 
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nent to constant density steady-states) at its intersection 
with the separatrix F. 

It is important to make, at this point, yet some ad¬ 
ditional remarks. First, the interval correspond¬ 

ing to the region CF in the bifurcation diagrams, was 
found to be unstable under variations of ur > 0. Corre¬ 
sponding (in)nstability regimes are illustrated in the top 
panel of Fig. where a Diagnostic Il-type diagram is 
shown, namely the bifurcation curve or — ||u(T'max)||a 
[continuous (red) line]. This diagram is plotted for fixed 
^max = 3000 and /? = 0.52 (recall that, in the previous 
case, for fixed (Jr = 0.01, it was found that /? = 0.52 e 
~ [0.5,0.53], i.e., in the chaotic regime); the rest 
of parameters are as in Fig. It is observed that for 
relatively small values of the SRS coefficient, namely for 
(Tfl < 0.03 (cf. grey-shaded area, labeled by SR), chaotic 
behavior persists. On the other hand, above this thresh¬ 
old, i.e., for cjR > 0.03, chaotic structures are destroyed, 
and the system enters into the metastability regime (la¬ 
beled by RW in the diagram). The ultimate steady-state 
is E3 for these parameters. Note that the instability of 
quasi-periodic and space-time periodic regimes under the 
influence of small increments of cr^j, occurs in a very sim¬ 
ilar manner, and can be plotted in similar bifurcation 
diagrams (results not shown here). 

Second, the interval Ip^K persists even in the absence of 
the rest of the higher-order effects, i.e., iov aR = n = v = 
0. This highlights the fact that the third-order dispersion 
plays a dominant role in the emergence of complex dy¬ 
namics. An example of the chaotic behavior, for /3 = 0.53 
and n = V = ur — 0, is shown in the bottom panels of 
Fig.0 In particular, the bottom left panel shows a part 
of a chaotic orbit for t G [600,650], of the 3D-projection 
of the flow on P3, for Xi = 5, X2 = 10, and X3 = 15. Fur¬ 
thermore, the three snapshots in the bottom right panel, 
show profiles of the solution corresponding to points PI, 
P2, and P3 of the chaotic path shown on the left, for 
t = 600, t = 625, and t = 650, respectively. 


IV. CONCLUSIONS 


namical picture is shared by various non-integrable per¬ 
turbations of Hamiltonian partial differential equations 
(such as the NTS, Sine-Gordon, and others), as these per¬ 
turbations may break the homoclinic structure of their 
integrable counterparts. However, the path to all the 
above dynamical scenarios can be traced in drastically 
different ways and to essentially distinct roots, even if 
the systems have similar origins for their dissipative na¬ 
ture manifested by the existence of an attractor, e.g., 
due to the presence of gain/loss, as in the case of CGL 
models. 

In particular, in our higher-order CGL model, keeping 
gain/loss - as well as other coefficients of the higher- 
order effects - fixed, we have shown that the competition 
between third-order dispersion and the SRS effect (in 
the presence of nonlinearity, dispersion, and gain/loss) 
can trace a path from Poincare-Bendixson-type behavior 
to quasi-periodic or chaotic dynamics. These dynamical 
transitions are also reminiscent of ones observed in fiber 
ring lasers, or in the path towards optical turbulence phe¬ 
nomena uniiszi- A conspicuous finding was that third- 
order dispersion chiefly appears to be playing a critical 
role in controlling the transition from periodic to quasi- 
periodic, and eventually to chaotic behavior, even in the 
absence of the rest of the higher-order effects. 

Our results highlight that higher-order effects may 
have a primary role for the birth of spatiotemporal tran¬ 
sitions in mixed gain/loss systems, suggesting further in¬ 
vestigations. First of all, in the framework of the model 
we considered herein, it would be particularly interest¬ 
ing to investigate more broadly the full six-parameter 
space, rather than its low-dimensional projection consid¬ 
ered herein. Furthermore, another interesting direction 
would be the identification of a low-dimensional attrac¬ 
tor, its dimension and dependence on the spatial length 
m, as well as the construction of the appropriate finite¬ 
dimensional reduced systems able to capture the effective 
low dimensional dynamics [38] . Lastly, it would also be 
interesting to investigate the role of higher-order effects 
in other autonomous systems with gain and loss. 


In conclusion, we have studied a physically impor¬ 
tant and broadly relevant higher-order Ginzburg-Landau 
equation, with zero diffusion. The considered model, is 
motivated by a higher-order nonlinear Schrodinger equa¬ 
tion, which finds applications in a variety of contexts, 
ranging from nonlinear fiber optics and deep water waves; 
the model also incorporates linear loss and nonlinear 
gain, while it is supplemented with periodic boundary 
conditions, which are relevant to optical cavities settings, 
such as ones employed, e.g., in ring lasers. 

Our analysis revealed that the infinite-dimensional dy¬ 
namics of this model can be reduced to a sequence of 
low-dimensional dynamical scenarios (fixed points, peri¬ 
odic and quasi-periodic, as well as chaotic orbits) that 
can be suitably revealed in reduced (two- and three- 
dimensional) phase space representations. Such a dy- 


Appendix A: Existence of a limit set (attractor) 


In this Appendix, we define an extended dynamical 
system associated to the initial-boundary value problem 
(©-Q- In particular, we briefly sketch the proof for the 
existence of a limit set-attractor, capturing all bounded 
orbits of this dynamical system, which initiate from suf¬ 
ficiently smooth initial data Q. 

The starting point of our proof is the power balance 
equation [55] : 


dt 


/ L pL pL 

jupdx = 27 / \u\'^dx + ‘2,5 / \u\‘^dx, (Al) 

-L J-L J-L 


satisfied by any local solution u G C([0,T], 

which initiates from sufficiently smooth initial data ug G 
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for fixed fc > 3. Here, denotes the 

Sobolev spaces of periodic functions [M] , in the fun¬ 

damental interval H = [—L,L]: 

= {u : H —> C, uand S j = 1,..., fc; 

u and dluior j = 1, ...,k — 1 , are 2 L-periodic}. 


Analysis of (Al I, results in the asymptotic estimate: 

limsup— [ \u{x,t)fdx<—^, (A2) 

t-^oo 2L J_j^ 0 


hence local in time solutions u € C{[0,T], are 

uniformly bounded in This allows for the defini¬ 

tion of the extended dynamical system 

ip{t, uo) : ip{t, uo) = u, 

whose orbits are bounded Vt > 0. Moreover, from the 
above asymptotic estimate, we derive, that if is 

endowed with the equivalent averaged norm 



then its ball 

B^{0,p) = {ueL^in) : \\u\\l<p^, 

attracts all bounded sets B G That is, there 

exists T* > 0, such that (p{t,B) C Ba, for all t > T*. 
Thus, we may define for any bounded set B G 
k > 3, its w-limit set in Lp‘{Vl)^ 


(where A is a real constant), which is an exact analytical 
solution of Eq. ([^ (for a MI analysis for the cw solution 
of Eq. 0 cf. Ref. [35]). This solution exists when the 
following dispersion relation holds: 

Wo = - kl/2 - pA^ko -I- i (7 -I- SA^) - A^, 

while A^ = — 7 / 5 , to suppress any exponential growth. 
This amplitude value is consistent with the equilibria 
(steady states) of the system. 

Now consider a small perturbation to this cw solution 

u{x, t) = [A -I- Ui{x, 

inserted into Eq. Q . Linearizing the system with respect 
to ui we obtain 


i(uit - kouix) - -uixx -I- A^('Ui -I- u*) = i(5A^(ui -I- u*) 

-I- i/3(3A:QUia; - 3ifcoMi XX "^Ixxx^ 

— ipA'^{ikoUi + ikoul + 2uix + 

- i[v - itTi^)A^(Mi 2 , -I- u\x), 

where star denotes complex conjugate. Solutions of the 
above equations are sought in the form: 

ui[xA) = 

where ci _2 are real constants, while k and w are the 
wavenumber and frequency of the perturbations. This 
way, we obtain the dispersion relation: 


^(s)=n 

s>0 

The closures are taken with respect to the weak topology 
of L^(H). Then, the standard (embedding) properties of 
Sobolev spaces imply that the attractor w(H) is at least 
weakly compact in L^(H), or relatively compact in the 
dual space For any initial condition 0, mq € 'Bj 

we denote its limit set by u}{uo) C oj{B). 

Appendix B: Modulational instability 


+ Pi{k)uj + p2{k) = 0 (B2) 

where 

Pi{k) = -2j3k^ + 2\-3fikl + ko + A^{2p + v - ifTfl)]fc, 
P2{k) = /3^fc® 

-I- [-3/3^A;g -I- /3/co - 2/3A'^{2p + ly - ian) - 1/4] 

-I- [9/3^/cq - 6/3^0 -I- A:q(1 - 6/3A^(^ + v - wr)) 

+ koA^{f3(6 — 6i5) + 3^ + 2^ — 2iaR) 

+ A^(ii5 -I- pA?{3p + 2v — 2iaR) — IjA:^, 


In this Appendix we provide the modulational insta¬ 
bility analysis of the cw state: 

u = u{t) = 9{x,t) = kox — ojot, (Bl) 


and it should be recalled that A^ = —y/il. It is clear that 
the system will always be modulationally unstable, since 
the solutions of Eq. (B2) are in general complex. 
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